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ABSTRACT 

In this paper we combine WFC3/UVIS F275W, F336W, and F438W data from the “UV 
Legacy Survey of Galactic Globular Clusters: Shedding Light on Their Populations and For¬ 
mation” (GO 13297) HST Treasury program with F606W, F625W, F658N, and F814W ACS 
archive data for a multi-wavelength study of the globular cluster NGC 6352. In the color- 
magnitude and two-color diagrams obtained with appropriate combination of the photometry 
in the different bands we separate two distinct stellar populations and trace them from the 
main sequence to the subgiant, red giant, horizontal and asymptotic giant branches. 

We infer that the two populations differ in He by AT = 0.029+0.006. With a new method, 
we also estimate the age difference between the two sequences. Assuming no difference in 
[Fe/H] and [a/ Fe], and the uncertainties on AT, we found a difference in age between the two 
populations of 10 ± 120 Myr. If we assume [Fe/H] and | ar/Fe] differences of 0.02 dex (well 
within the uncertainties of spectroscopic measurements), the total uncertainty in the relative 
age rises to ~300 Myr. 

Key words: stars: Population II — globular clusters individual: NGC 6352 


1 INTRODUCTION 

The Hubble Space Telescope ( HST) “UV Legacy Survey of Galac¬ 
tic Globular Clusters: Shedding Light on Their Populations and 
Formation” (GO 13297, PI Piotto) is a Treasury program to study 
globular cluster (GC) stellar populations using multi-wavelength 
high-precision HST photometry and astrometry. The program is 
fully described in IPiotto et alj 12015L hereafter Paper I). The main 
purpose of this survey is to identify multiple stellar populations 
and study their properties, including their relative ages, chem- 


* Based on observations with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS 5-26555 under GO 13297. 
t E-mail: domenico.nardiello@studenti.unipd.it. 


ical compositions, spatial distributions, and kinematics. Within 
GO 13297 we collected F275W, F336W, and F438W WFC3/UVIS 
images, approximately overlapping the F606W and F814W data 
from GO 10755 (PI Sarajedini, see lsaraiedini et al]l2()07l) . 


In this paper we analyze photometry of stars in the Galactic 
GC NGC 6352. Our main purpose is to infer the helium content 
and determine the relative ages of its stellar populations. NGC 6352 
is a Bulge (Galactic coordinates: l = 341, b = -7.2), metal- 
rich GC, located at 5.6 Kpc from the Sun [(m - Af)y=14.43, 
E(B- V)=0. 22.lHarrisl 19961 2010 revision], and having luminosity 
M v = -6.47 l Harris! 19961 2010 revision) and mass of 3.7 x 10 4 M,., 
l iMarks & Kroupal l201dl) . From high-resoluti on spe ctr oscopy of 
nine horizontal branch (HB) stars, iFeltzing, Primas & Johnsonl 
d2009t) confirmed the results of previous works, that this clus- 
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ter has a high metallicity ([Fe/H]~ -0.55) and is a-enhanced 
([a/Fe]~0.2). 

The paper is organized as follows. Section 2 presents the data 
and data analysis. Section 3 and 4 show the characteristics of the 
two stellar populations hosted by NGC 6352, and their properties 
as seen with different combinations of colors and magnitudes. In 
Section 5 the difference in the helium content of the two stellar 
populations is calculated. In Section 6 we describe the new method 
we developed to estimate the difference in age between the two 
stellar populations. In Section 7 there are the conclusions. 


2 OBSERVATIONS AND DATA ANALYSIS 

In order to identify the multiple stellar populations in NGC 6352 
we used WFC3/UVIS images of GO 13297 and ACS/WFC data 
of GO 10775. The WFC3 data set consists in 2x706 s +2 x 800 s 
F275W, 4 x 311 s F336W, and a 58 s + 72 s F438W images. ACS 
data include 4 x 140s + 7s F606W and 4 x 150s + 7s F814W 
images, overlapped to WFC3 images. We also reduced ACS im¬ 
ages collected in F625W and F658N bands within GO 12746 (PI 
Kong). Exposure times are 2x 150 s in F625W and a 650 s + 643 s 
in F658N. A detailed description of the reduction of GO 13297 
data is provided in Paper I (see their Sect. 5). We ha ve used t he 
photom etric and astrometric catalogs published bv lAnderson et alj 
( 120081) for GO 10775 data, while photometry and astrometry of 
archive ACS /W FC i mages from GO 12746 have been carried out 
as in lAnderson et alj l l2008!) . 


3 THE COLOR-MAGNITUDE DIAGRAMS OF NGC 6352 

In Paper I (see their Fig. 2), we have shown that the color- 
magnitude diagram (CMD) of NGC 6352 is not consistent with a 
single stellar population. The mp 8 i 4 w vs. m F 606 w - »J F si 4 w CMD of 
the cluster members of NGC 6352 is shown in Fig.Q] in black there 
are the stars that, on the basis of their proper motions, have high 
probability to be cluster members, in light blue the rejected stars. 
To derive stellar proper motions, we determined the displacement 
between the stellar positions in GO 10775 data-set (2006.4) and in 
GO 13297 data-set (2013.7 and 2014.5) by following the method 
described in detail by [Piotto et aD d2012L see their Sect. 4). The 
maximum time baseline for the proper motion measurements is 8.1 
yrs. 

Figure Eh shows the '«F606 W vs. 
CF275W,F336W,F438W=( m F275W _m F336w) _ ( m F336W _ ' w F438w) di a_ 

gram for proper-motion-selected cluster members. Asymptotic 
Giant Branch (AGB), stars selected from the rap 336 W vs. ;«f 336 w- 
' JI F8i4w CMD. are represented with starred symbols. Two distinct 
sequences are clearly visible and can be continuously followed 
from the main sequence (MS), the sub-giant branch (SGB), the red 
giant branch (RGB), to the HB, and the AGB. Panels b-e of Fig. [2] 
show a collection of CMDs for NGC 6352 and reveal the pattern 
of its multiple sequences. 


3.1 Multiple stellar populations in NGC 6352 

In Fig.[3]we use the m F3 36w- m F438w vs. m F2 75w- m F336w two-color di¬ 
agram to separate the two populations, named population-a (POPa) 
and population-b (POPb) hereafter. Panel (a) of Fig. [3] shows the 



Figure 1 . The /« F 8i4W vs. m F 606W — ^fsmw diagram of NGC 6352 . In black 
and in light blue are plotted cluster members and field stars, respectively. 
Cluster members have been selected on the basis of their proper motions. 
The vector-points diagram of proper motions is plotted in the inset. Black 
points within red circle and light blue points outside red circle indicate 
members and field stars, respectively. 


m F 8 i 4 w vs. m F 606 w-ntF 8 i 4 w CMD, corrected for differential redden¬ 
ing and only for cluster-member stars, selected using proper mo¬ 
tions. The green dotted lines identify the four regions in the CMD 
which include MS, SGB, RGB and HB stars. Again, AGB stars are 
plotted using starred symbols. All stars that, on the basis of their 
position on the CMD, are possible binaries, blue stragglers or sur¬ 
vived field stars have been excluded from the following analysis 
and plotted as gray crosses. 

Our recent papers demonstrated how two-color and color- 
magnitude diagrams made with appropriate combination of far- 
ultraviolet, ultraviolet and blue magnitudes represent a very ef- 
ficient tool to iden ti fy mu ltiple stellar populations in a GC 
dMilone et al J|20 1 2j] |cl l20 1 3l : PaperI). Panels (bl), (b2), (b3), and 
(b4) of Fig.[3]show the /K F 3 36W -ra F4 3 8W vs. m F 275 w m F 336 w two-color 
diagrams for the four evolutionary sequences highlighted in the 
left-panel CMD. Two sequences are clearly visible in each dia¬ 
gram. We drew by hand a straight continuous line to separate the 
two groups of POPa and POPb stars in the SGB, RGB + AGB, 
and HB, and colored them green and magenta, respectively. The 
two sequences of AGBa and AGBb stars are separated by the black 
dashed line in panel (b2). In order to separate the two MSs we used 
the m F 8 i 4 w vs. C F 275w,F336w,F438w diagram plotted in panel (cl) of 
Fig. [3] Indeed, in this pseudo-CMD the double MS of NGC 6352 is 
better distinguishable than in the m F 336 w ,M F 438 w vs. m F 275 w _m F 336 w 
plane. 
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Figure 2. Overview of the main features in the CMDs of NGC 6352. Panel (a) shows the mp 606 W vs. Cf275W,F336W,F438W pseudo-CMD, in which the split of 
all the sequences is clear (starred points represent the AGB stars); panel (b) shows the region around the HB in the mp 275 w vs. m F336W ~ m F438W CMD; panel 
(c) highlights the SGB split in the mp336W vs. mp275w _ '”F336W diagram; panel (d) shows the RGB region in the mp8i4W vs. Cf336W,F438W,F814W pseudo-CMD. 
The Hess diagram of the MS stars shown in panel (a) is plotted in panel (e) and highlights the two distinct sequences. 


The green and the magenta fiducial lines, superimposed on 
the MS, are the fiducials of the MSa and the MSb, respectively, 
and have been obtained by using the following procedure. As a 
first guess, we have selected by eye a sample of MSa and MSb 
stars and derived for each of them a fiducial line by fitting a spline 
through the median values of Cf275w,f336w,f438w obtained in suc¬ 
cessive short intervals of magnitude. W e iterate d this step with 
a sigma-clipping procedure. Then, as in iMilone et alJ 00151 . Pa¬ 
per II) , we verticalized the MS in such a way that the green 
and magenta fiducials translate into vertical lines with abscissa 
=0 and -1, respectively. The histogram of the dis- 
tribution in A£ plotted in panel (c2) is clearly bimodal. 

CF275W.F336W.F438W J 


We considered stars with A? > -0.5 as MSa members 

CF775W.F336W.F438W 

and the remaining MS stars as MSb objects. Panel (c3) shows the 
"JF 814 W vs. A5 verticalized diagram for the MS stars. 

We colored MSa and MSb stars in green and magenta, respectively. 
These colors will be consistently used hereafter. 

3.2 The chemical composition of the HB stars 

Typically, multiple sequences along the CMD correspond to dis¬ 
tinct stellar populati ons with different content of helium and light 
elements ( see, e.g., Marino et a O 120081; 1 Yong & Grundahil 120081 : 
Paper II). iFeltzing, Primas & Johnsonl J2009I) measured chemical 
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Figure 3. Procedure for the selection of POPa and POPb stars. Panel (a) shows the mpsnw vs. hif 606W - »tF8l4W CMD of NGC6352: gray crosses are the 
stars that are possible binaries and blue stragglers and that are not used in this work. Panels (b) show the ;«F 336 W — »>F 438 W vs. mp 275 W - m / ; 336 tv diagrams for 
MS, S GB. R GB, AG B and HB stars of POPa (green) and POPb (magenta). In panel (bl) we also plot the stars in common with the spectroscopic catalog of 
IFeltzing. Primas & lohnsorj 1200 91). Panels (c) show the procedure used for the selection of MSa and MSb stars (see text for details). 


abundances for nine HB stars of NGC 6352 from high signal-to- 
noise UVES@VLT spectra. They have confirmed that NGC 6352 
is a metal-rich GC ([Fe/H]=-0.55±0.03), and is enhanced in a 
elements by [ar/Fe]~0.2 dex. Feltzing and collaborators also de¬ 
tected significant star-to-star sodium variation in close analog y with 
what is observed in most M ilky-Way GCs (e.g., Kraft et ali' 19931 : 
iGratton, Sneden & Carrettal2004) . 


In order to investigate the chemical content of POPa 
and POPb , we have _exploited the spectroscopic results by 
IFeltzing, Primas & Johnson ( b009l) . Our photometric catalog in¬ 
cludes seven stars studied by Feltzing, Primas & Johnsoij l l2009l) . 
Five of them belong to the HBa and are Na-poor ([Na/H]< -0.38). 
The other two are Na-rich ([Na/H]> -0.27) and are HBb members. 
The two groups of Na-rich and Na-poor stars both have the same 


mean metallicity [Fe/H]= -0.55 ± 0.02. In panel (bl) of Fig. [3] 
we plotted these HB stars as triang les w ith the corre sponding ID 
number adopted bv lFeltzing, Primas & Johnson! J2009I ): orange tri¬ 
angles are for Na-poor stars, cyan triangles show the Na-rich ones. 


4 MULTI-WAVELENGTH VIEW OF MULTIPLE 
POPULATIONS 

In the previous section we used the m F si4w vs. Cf275w,f336w,f438w 
pseudo-CMD and the otf336w-'Uf438w vs. t«f275w-'«f336w two-color 
diagram, where multiple sequences are clearly distinguishable, in 
order to select the members of the two stellar populations along 
the MS, the SGB, the RGB, and the HB of NGC6352. We can 










































now combine all different bands in order to study the behavior (e.g. 
color differences) of the two populations in all possible CMDs we 
can construct with our data. An example is shown in Fig. [4] Note 
that POPa stars are redder than POPb ones in the m F2 75 W vs. m F 275 w- 
'UF 336 W CMD, while they are bluerjn the < n F 336 w vs. m F 336 w-'UF 438 w 
CMD, as in the case of 47 Tuc l lMilone et alj|2012<3) . The two up¬ 
per panels of Fig.[5]show the m F si 4 w vs. m x - »Jf 814 w fiducial lines 
for the RGBs (top panels) and MSs (middle panels) of the two pop¬ 
ulations, where X=F275W, F336W, F438W, F606W, F625W, and 
F658N. 

POPa is redder than POPb in most of the CMD s of Fig. [5] 
in close analogy with what is obser ved in most GCs d Be 11 i n i et alj 
120101 : iMilone et alj l2012al lcl 1 201 3l ). In visual and ultraviolet the 
m x - m F 8 i 4 w color separation between POPa and POPb increases 
with the color baseline, and is maximum for m F275W ~ m F8i4w- The 
fiducial of the two populations are almost coincident in m F336W - 
'UF 814 W, where POPa is slightly bluer than POPb. 

The procedure used to dete rmine each fiducial line is based on 
the naive estimator dSilvermanll 19861) . First, we defined a series of 
N mpsi 4 w magnitude intervals with a given width (w). Magnitude 
intervals are defined over a grid of N points separated by step of 
fixed magnitude (s). We calculated the median magnitude m F8 i 4 w,i 
and color (m x - m F 8 i 4 w)j of stars within each magnitude interval 
i = 1 These median points have been smoothed using the 

smoothing technique of boxcar averaging, in which each point of 
a vector is replaced with the average of the M adjacent points. Fi¬ 
nally, the smoothed median points were interpolated with a cubic 
spline. 

The choice of w is the result of a compromise. On one hand, 
we require small magnitude intervals to account for the detailed 
structure of the fiducial line. On the other hand, we need a large 
width to include in the magnitude interval a large number of stars 
for a statistically significant measurement of the median color and 
magnitude. We used different values of w for stars with different 
luminosity. Specifically, we assumed w = 0.5 for RGB stars with 
15.00 < m F 8 i 4 w < 17.00, w = 0.25 for SGB stars with 17.00 < 
m F 8 i 4 w < 17.75, and w = 0.05 for MS stars with 17.75 < m F si 4 w < 
19.00. In all the cases we used s = w/ 3. 

In order to estimate the error associated to the colors (<Xfjdcoi) 
and magnitudes (cr fldmag ) of each observed fiducial we used the fol¬ 
lowing procedure. For each stellar population selected in Fig. [3] 
we assigned a subsample of artificial stars (A Ss) and randomly ex - 
tracted a star from it by following the recipe in lMilone et alj j2009h . 
The AS subsample consists of all the ASs with similar magnitudes 
(within 0.3 m F gi 4 w magnitudes) and radial distances (less than 100 
pixels from the observed star). This method produces a catalog of 
simulated stars with almost the same luminosity and radial distri¬ 
bution of the observed catalog. We applied the procedure described 
above to estimate the fiducial line of the sample of AS stars, and 
calculated the difference between the fiducial of the ASs and the 
real stars. We repeated this step 100 times. The color error asso¬ 
ciated to each point of the observed fiducial is calculated as the 
68.27 th percentile of this distribution. 


5 THE HELIUM ABUNDANCE OF THE TWO STELLAR 
POPULATIONS 

The helium abundance of a stellar population is a fundamental 
ingredient to understand its evolution, and to shed light on the 
chemical-enrichment and the star-formation history in a GC. In ad¬ 
dition, accurate helium determination is crucial to estimate several 
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parameters of the host star cluster, like age, mass, and mass func¬ 
tion. 


is only feasible in rare cases and for a tiny subset of stars (e.g., Behi 

2003 

; Villanova. Piotto & Gratton 2009 

Dupree, Strader & Smith 

2011 

: Pasquini et al. 2011; Marino et al 

2013). On the other hand. 


the method based on multi-wavelength photometry of multiple se¬ 
quences that we developed can be applied to all clusters and pro¬ 
vides a reliable estimate of the He abundance differences among 
the different populations and c an reach int ernal e rrors smaller 
then spectroscopic methods (see IMilone et al J[2012d : Paper II and 
references therein). In the few cases where helium has been in¬ 
ferred from both spectroscopy of HB stars and photometry, the 
results from the t wo techn i ques are in fairly agreement. In the 
case of NGC 2808, iMarino et alj J20l4 i) derived an helium abun¬ 
dance T = 0.34 ± 0.01 ± 0.05 (internal plus systematic uncer¬ 
tainty) from spectroscopy of HB stars slightly bluer than the RR 
Lyrae instability strip. This result is consistent with the helium 
abundance inferred from the analysis of multiple MSs, where 
the middle MS (which would be associated to the stars an alyzed 
by Marino and coll aborators), has Y ~ 0.32 JPiotto et alj|2007l : 
iMilone e t all201 2bl). S imilar ly in the case of N GC 6397, both pho- 
tometr y ddi Criscienzo. D’ Antona&A^ntura J20ld: IMilone et alj 
l2012al) and spectroscopy i Mucciarelli et al. 20141) conclude that 
second-population stars ar e slightly e nha nced in Y by AT ~ 0.01 
dex. In the case of M 4, Ivillanova et alj J 20 T 2 I) found that the sec¬ 
ond population is enhanced in helium by 0.04-0.05 dex with re¬ 
spect to the first one . This resul t is in apparent disagreement with 
results bv lNardiello et alj d2015h who inferred an helium difference 
of 0.020 ± 0.008 dex between the two MSs of this clusters. How¬ 
ever, it should be noted that the discrepancy could be due to NLTE 
that affect the helium line at 5875.6 A (see e.g. Marino et al. 2014, 
Sect. 4 for the case of NGC 2808). 

Helium has been also inferred in RGB stars from near-infrared 
chromospheric transition of He I at 10 830A in NGC 2808 and 
to Centauri. The spectra suggest helium abundances of T <0.22 
and T =0.39-0.44 ( AT > 0 .17) for the two analyzed stars in 
to Centauri jDupree & Avrettll2013l). while Na-rich and Na-poor 
stars of NGC 2808 analyzed by lPasquini et al J d201ll) differ in he¬ 
lium by AT >0.17, with the sodium-rich star being also helium 
rich. Such helium variation are larger than those obtained from the 
analysis of multiple sequences in the CMD, that a r e ~ 0.07 for 
NGC 2808 dD’Antona et alJl2005l :JPiotto^t_alj|2007l: IMilone et al] 
l2012bl) and ~ 0.13 for o> Centauri dKing et al.ll2012l) . 

Figure [5] illustrates the procedure to infer the helium content 
of the stellar populations in NGC 6352. Upper panels show a col¬ 
lection of m F8 i 4 w vs. m x - m F8 ] 4W fiducial lines, where X=F275W, 
F336W, F438W, F606W, F625W and F658N. We have used green 
and magenta continuous lines to plot the fiducials along the RGB 
(top panels) and the MS (middle panels) of the two populations. We 
have calculated the color difference between the green and the ma¬ 
genta fiducials at the reference magnitudes indicated by the black 
dashed lines at m F8 J 1 4 W =18.75 for the MS, and 16.50 for 

the RGB. 

Results are shown in the lower panels of Fig. [5] We found that 
A (m x - m F 8 i 4 w) has negative values only for X=F336W. When the 
other filters are used, the m x - m F8 i 4 w color separation between 
the two RGBs and MSs is maximum for X=w F2 75 W and decreases 
towards redder wavelengths of the X filter. 

In order to estimate the effective temperatures (T cff ) and grav¬ 
ities (logg) of MS and RGB stars with «J F8 i 4 w =m F8 i 4 W we used the 
isochrones bv lDotter et alj J2008I) that best match the CMD. 
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Figure 4. The mp 275 W vs. »iF 275 W-«iF 336 W (left panel) and ;«f 336 W vs. "iF336W-”*F438W (right panel) CMDs. In the first CMD, POPa stars (green dots) are 
redder than POPb stars (magenta crosses), while in the /Wf 336 W vs. »iF336W _m F438W CMD POPa stars have, on average, bluer colors than POPb stars. 


— 

m F 814 W 

AmF 275 W,F 814 W 

A/nF 336 W,F 814 W 

A/HF 438 W.F 814 W 

A»fF 606 W,F 814 W 

An!F 625 W,F 814 W 

A»lF 658 N,F 814 W 

19.50 

0.300±0.028 

-0.016±0.017 

0.049±0.011 

0.012±0.003 

0.005±0.004 

0.001±0.004 

19.25 

0.234±0.026 

-0.009±0.012 

0.045±0.007 

0 .011±0.002 

0.006±0.003 

0.005±0.003 

19.00 

0.204±0.019 

-0.003±0.009 

0.048±0.006 

0 .012±0.002 

0.009±0.002 

0.005±0.003 

18.75 

0.174±0.014 

-0.007±0.007 

0.051±0.005 

0 .010±0.002 

0.009±0.002 

0.004±0.002 

18.50 

0.128±0.011 

-0.011±0.006 

0.045±0.004 

0.007±0.002 

0.005±0.002 

0 .001±0.002 

18.25 

0.111±0.008 

-0.008±0.004 

0.038±0.003 

0.008±0.001 

0.005±0.001 

0 .001±0.002 

17.00 

0.168±0.038 

-0.056±0.030 

0.052±0.018 

0.016±0.008 

0.016±0.008 

0.022±0.009 

16.50 

0.179±0.028 

-0.039±0.016 

0.046±0.008 

0.013±0.005 

0.008±0.005 

0.009±0.005 

16.00 

0.158±0.030 

-0.033±0.020 

0.044±0.008 

0.011±0.005 

0.010±0.007 

0.010±0.006 


Table 1. Color difference Amx,F8l4W at different 


Note that we are mainly interested in relative ages rather than 
in absolute ones. Indeed, as demonstrated in next section, accu¬ 
rate relative ages are feasible using our dataset. Here, we just need 
reference models to calculate the color differences of the two pop¬ 
ulations. As discussed in Section [6] the best match with observa¬ 
tions could be obtained using isochrones with [Fe/H]=-0.67, in 
fair agreem ent with the value inferred from high-resolution spec¬ 
troscopy by ICarretta & Grattonl d 1997lf . [a/Fe]=0.4, and an age of 
13.1 Gyr [adopting (m - M) v = 14.43, an E(B - V) = 0.22. IHarrisl 
119961 2010 edition]. We assumed Y = 0.256 for POPa, and dif¬ 
ferent helium abundances with Y ranging from 0.256 to 0.330, in 
steps of 0.001, for POPb. Since there are no measurements of C, N, 
O abundance for NGC 6352 stars, we assumed for POPa and POPb 
the same C, N, O content as the two main stellar populations of 
47 Tuc. Specifically, we used for POPa [C/Fe]=-0.20, [N/Fe]=0.20, 
[O/Fe]=0.40, and assumed that POPb stars are depleted in carbon 


and oxygen by 0.20 and 0.15 dex, respectively and nitrogen en¬ 
hanced by 0.7 dex with respect to POPa stars. 

We used the ATLAS 12 code l lKurucziri993l ; ISbordone et ai] 
l2004ll to calculate atmospheric models for the MS and RGB stars 
with mFSHW^m S^w, while synthetic s pectra have been generated 
with S YNTHE dKurucz & Avrettll98l!) with a resolution of R=600 
from 2000 to 10000A. 

Synthetic spectra of POPa and POPb stars have been inte¬ 
grated over the transmission curves of the WFC3/UVIS and the 
ACS/WFC filters used in this paper in order to determine synthetic 
colors. The colors from synthetic spectra with different Y have been 
compared with observations in order to estimate the helium abun¬ 
dance of POPb. 

The best fit between observed and theoretical colors at 
m F8i4w =18-75 can be obtained assuming POPb enhanced in he¬ 
lium by AY=0.033 dex. We have repeated this procedure for 
m F8 U i4w =18-25, 18.50, 19.00, 19.25, and 19.50 for the MS, and 
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Figure 5. Top Panels: RGB fiducial lines for POPa (green solid lines) and POPb (magenta dotted lines) in six H2F814W vs. m\ - JHF814W CMDs (where 
X=F275W, F336W, F438W, F606W, F625W, F658N); the panel insets show the color distance between the two fiducials at 'Wpsutw = 16.5. Middle Panels: 
same as top panels, but for the MS; the color difference indicated in the inset is measured at 'AAAy = 18.75. Bottom panels: m\ — mp 8 i 4 w color distance 
between MSa and MSb at nipsilw = 18.75 (left) and between RGBa and RGBb at HipsMW = 16.50 (right) as a function of the central wavelength of the X 
filter. 


m F 8 i 4 w =16.00, 16.50 and 17.00 for the RGB. The results are listed 
in TableUI The values of r c jf, log(g), and AY that provide the best fit 
to the observed color differences are listed in Table[2]for each value 
of »Jf 814 w- avera g e > we have AY=0.029±0.006, where the er¬ 
ror represent the 68.27th percentile of the distribution of the sorted 
residuals from the mean value. To estimate AY we only used visual 
bands (F606W, F625W, F658N, and F814W), because they are not 
affected by light-elements variations. For this reason, the final re¬ 


sult is not conditioned by the assumptions on the C, N, and O abun¬ 
dances. In the case where C, N, and O abundances were known, the 
final value of AY, obtained using all the photometric bands, would 
not change. As a test, we repeated the procedure described above 
by assuming different values of C, N, and O: we found the same 
result. 
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m CUT 

"T 814 W 

TEFF,POPa 

TEFF,POPb 

log gPOPa 

log gPOPb 

AT 

19.50 

3673 

3677 

2.51 

2.53 

0.018 

19.25 

3678 

3684 

2.64 

2.69 

0.023 

19.00 

3687 

3692 

2.86 

2.91 

0.030 

18.75 

3693 

3698 

3.06 

3.10 

0.033 

18.50 

3698 

3702 

3.23 

3.26 

0.024 

18.25 

3700 

3705 

3.37 

3.40 

0.031 

17.00 

3766 

3767 

4.25 

4.26 

0.042 

16.50 

3764 

3766 

4.29 

4.30 

0.026 

16.00 

3754 

3756 

4.40 

4.41 

0.030 

0.029±0.006 


Table 2. Parameters used to simulate synthetic spectra of POPa and POPb 
stars and estimation of helium difference between the two populations for 
different m™^ w . 


6 RELATIVE AGES OF THE TWO STELLAR 
POPULATIONS 

The relative ages of the multiple stellar populations hosted by GCs 
is an important issue to understand how GCs formed. We can spec¬ 
ulate that the phenomenon of multiple stellar populations in GCs is 
due to the presence of different generations of stars, formed in dif¬ 
ferent epochs: a first stellar generation (FG) characterized by pri¬ 
mordial helium and chemical composition similar to that of field 
stars with the same metallicity, and a helium-enhanced second gen¬ 
eration (SG), characterized by stars depleted in C and O and en¬ 
hanced in N and Na, born from material processed at high tem¬ 
perature b y FG stars. If we exclude the “an omalous” GCs, suc h 
as oj Cen i Nom^& DaCostal 1 1995h . M22 l lMarino et alJl20lTh . 
and M2 I Yon si et al 2014 . Paper II), stars of different popula¬ 
tions in “normal” G Cs show negligible differences in metallicity 
dCarretta et alj|2009l ). Because of this, we can exclude supernova 
ejecta as possible polluters of the material from which SG stars 
formed. 

Current scenarios for the formation of multiple stellar pop- 
ulatio n s in GCs pred ict different timescales. The AGB scenario 
to’ Antona et all2002l ) predicts that envelopes of intermediate mass 
AGB stars are the cause of pollution. In models aimed at repro¬ 
ducing the observed chemical patterns in “normal” globular clus¬ 
ters with AGB ejecta, the second generation star formation epoch 
exten ds between about 40 and 90-100 Myr to’Ercole et alJl201~oL 
l20H) . Another suggested mechanism is the formation of stars from 
material ejected by Fast-Rotating Massi ve Stars (FRMS) during the 
phas e of core H-burning toecressin_etaljj2007h or massive bina¬ 
ries dBastian et alj 120131 : ISalaris & Cassisil 20141) . In these cases, 
the timescale for the formation of the seco nd gen eration must 
be of the order of a few Myrs or even less dBastian et alj 1 201 3l ; 
ISalaris & Cassisill2014l : ID’ Antona et al.|[2QI4l) . 

In this work we present a method to set upper limits on the dif¬ 
ference in age between the two populations of NGC 6352. For the 
first time, we estimate the relative age of two stellar populations in 
a normal GC, i.e. with no internal variation in [Fe/H]. Note that we 
are interested in measuring relative ages. The corresponding abso¬ 
lute ages should be only regarded as indicative (depending on the 
adopted distance modulus and reddening, as well as the reference 
model). 

We consider a set of isochrones characterized by [Fe/H]= 
-0.67 and [a-/Fe]=0.4, that are the values that best fit all the stars 
in the /mfsmw versus /mf 606 w - m F 8 i 4 w> as explained in Sect. [5] 

We also used the m F 8 i 4 w versus m F 606 W - m F 8 i 4 w CMD for 
computing the relative ages between the two populations, because 
in these filters the effect of light elements variation is not signifi¬ 


cant dSbordone et al1 l20 1 lb iMilone et alj|20l2cl) . The MS turn-off 
is the classical age indicator of a simple stellar population. For this 
reason, we considered a region of the CMD around the MS turn-off 
of the two populations [panel (a) of Fig. [6). We obtained fiducials 
and the associated errors in color and magnitudes (colored regions 
between dashed lines) for POPa (green) and POPb (magenta) using 
the same procedure described in Sect. [4] 

In order to compute the relative age between the two pop¬ 
ulations, we treated POPa and POPb as simple stellar popula¬ 
tions and measured their ages independently. The age of each 
population is estimated by comparing the observed m F 8 i 4 w versus 
«F 606 W - '«F 8 i 4 w CMD with a grid of synthetic CMDs with the 
same [Fe/H], [ff/Fe] and T, and ages that run from 12000 Myr to 
15000 Myr, in step of 50 Myr. 

For each isochrone a synthetic CMD was built in this way: 
we obtained a first guess synthetic CMD by interpolating the 
isochrone on the observed magnitudes, in such a way that to each 
' JI F 8 i 4 w magnitude in the catalog of real stars is associated a color 
' j; F606w -m F8i4w on the isochrone As second step, we broadened this 
synthetic CMD by adding to the color and magnitude of each syn¬ 
thetic star a random Gaussian error, with a dispersion equal to the 
error of the associated real star. The final result is a synthetic CMD, 
characterized by the same number of stars and a similar luminosity 
function to that of the observed CMD. 

Panels (b) and (c) of Fig. [6] show an example of isochrones 
and corresponding synthetic CMDs. For each synthetic CMD we 
computed the fiducial line using the same procedure as for the ob¬ 
served data. For each step in age, we built N = 100 synthetic simple 
stellar population CMDs, obtained the fiducial lines and computed 
the average of all the 100 fiducial lines in order to obtain the final 
synthetic fiducial line. 

We compared the observed fiducial line of each population 
with the synthetic ones. This method is more robust than the di¬ 
rect comparison between the observed CMD and the isochrones, 
because we compared observed and theoretical fiducial lines which 
were computed in the same way and are affected by the same sys¬ 
tematic errors introduced by the method (as, e.g., the error due to 
the smoothing). For each step in age, we calculated the x 2 of the 
best fit between the synthetic and the observed fiducial, in a mag¬ 
nitude interval 17.00 < otfshw < 18-2. We fitted thex 2 as a func¬ 
tion of the age with a second-order polynomial, in order to find 
the age that minimizes the^ 2 . This fit was performed in two steps: 
first, we fitted all the points with a parabola to obtain a first-guess 
value of the age, f 1 , that minimizes the^ 2 [panels (di) and (ei) of 
Fig- ED- Then we considered only the points inside an interval be¬ 
tween t 1 - 500 Myr and t l + 500 Myr and fitted a new second-order 
polynomial [panels (di) and (e 2 )] to find the best value that mini¬ 
mizes the^ 2 . 

We independently run this procedure for POPa and POPb. For 
POPa we used a set of isochrones calculated for [Fe/H]= -0.67, 
[a/Fe]= +0.4 and Y = 0.256, while for POPb we used isochrones 
with the same metallicity and or-enhancement, but enhanced in he¬ 
lium by AT = 0.029 (i.e. with T = 0.285) as obtained in Sec- 
tion[5] The best fit gives an age of (13580 ± 80) Myr for POPa and 
(13570 ± 80) Myr for POPb. We want to emphasize the fact that 
absolute ages are not significant in this context, and that we are in¬ 
terested only to relative ages. For this reason, hereafter we provide 

only relative ages. The age error is given by a = ^crp H0T + o"p IT , 
where <r PH0X is the error due to photometric errors and <t F it is the 
error on the fit of the^ 2 values. 

The error due to photometry has been computed by perform- 
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Figure 6. Procedure used to obtain the relative age of the two populations. In panel (a) the observed mp8i4w vs. mp606W ~ '”F814W CMD is shown. In green 
and magenta are the fiducial lines of the POPa and POPb, respectively. Panel (b) shows the best-fit isochrones for the two populations, while in panel (c) we 
plot the synthetic CMD and the fiducial lines built using the isochrones calculated for the two population. Panels (d) and (e) illustrate the procedure adopted 
to derive the age of POPa and POPb, respectively. In panels (di) and (ei), filled circles indicate the values ofy ' 2 for different ages while the dashed lines are 
the first-guess second-order polynomial fit. The insets ( 82 ) and fe) show the best-fit second-order polynomial for POPa (green) and POPb (magenta, see text 
for details). 


ing the same analysis as described above and shifting the observed 
fiducial line in color by a quantity equal to ±cr fldcol and in mag¬ 
nitude by a quantity equal to ±<x fldmag . To estimate <xfit, we used 
the bootstrap re-sampling of the data to generate 1000 samples ran¬ 
domly drawn from the original data sets; for each sample, we cal¬ 
culated the age that minimize they' 2 as above described. Finally we 
computed the mean age and its standard deviation, and adopted the 
latter as the uncertainty in the x 2 fit- In conclusion, the two pop¬ 
ulations have a difference in age Aage = AGEpop a - AGEpopb = 
(10 ± 110) Myr. 


Figure [6] shows the whole procedure: in panel (a) we plot the 
CMDs of POPa and POPb (light green and magenta) and their 
fiducial lines. In panel (b) the best-fit isochrones are plotted for 
each population; these are the isochrones for which the corre¬ 
sponding synthetic fiducials, shown in panel (c) with the synthetic 
CMDs, give the minimum x 2 - Panels (d) and (e) show the procedure 
adopted for the fit of theyf 2 for POPa and POPb respectively. In pan¬ 
els (di) and (eD are shown the first guess fits (blue dashed lines), 
while in panels (d 2 ) and (eg) show the final fits (green and magenta 
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curves for POPa and POPb respectively). The vertical dashed lines 
represent the value of the age that minimizes the^ 2 . 

We performed the same analysis using CMDs with different 
color bases! /uf 625 w — ^F 8 i 4 w and /Upg 5 gj^ — ^F 8 i 4 w* These CM13s are 
characterized by larger photometric errors. We found Aage = 80 ± 
140 Myr in the case of the mpsww vs. mp 625 w - /kfsmw CMD, and 
Aage = 220 ±210 Myr in the case of the mp 814W vs. m F6SSN - m F8 i 4W 
CMD. The results are consistent with what found using the ;7 Jfsi4w 
versus mF 606 w - WF 8 i 4 W> but errors due to photometric uncertainties 
are larger. For this reason we decided to carry out further analysis 
only using the mp 8 i4 W vs. m F60 6w - '«F8i4w CMD. 

We tested whether our results depend on the adopted magni¬ 
tude interval for the fit. We iterated the procedure described above 
by changing the starting point of the magnitude interval between 
m F 8 i 4 w = 17.80 and mp 8 i 4 w = 18.10, and the ending point be¬ 
tween /77F814W = 17.00 and mpsnw = 17.30, with steps of 0.05 mag. 
Within 50 Myr, we found no difference in the average age. Such age 
difference is consistent with zero within our measurement errors. 

We tested our results using a different dataset of isochrones, to 
prove that results are indep endent from adopted models. We con¬ 
sidered BaSTI isochrones dPietrinferni et alj|2004 l2009h and we 
performed the same analysis described above. We obtained a dif¬ 
ference in age between the two populations of (-80 ±110) Myr, in 
agreement with the results obtained previously. 

There is some dispersion among the abundance measurem ents 
of N GC6352 in the literature (see iFeltzing. Primas & Johnson! 
120091 and discussion therein). Therefore, we re-iterated the 
same analysis using another set of isochrones with differ¬ 
ent [Fe/H] an d Itt/Fel. Following the spectroscopic results 
by IFeltzing, Primas & Johnsonl J 2 OO 9 I) . we considered a set of 
isochrones with [Fe/H]=-0.55 and [o-/Fe]=+0.2 and performed the 
same procedure as described above. We found for the two pop¬ 
ulations a slightly different relative age (Aage = 70 ± 110 Myr). 
We thus confirm that the two populations have a difference in age 
within ~ 110 Myr and that the adopted values of [Fe/H] and [or/Fe] 
do not change our conclusions. 

In addition to the photometric error, we also considered the 
uncertainty in AT, that in the previous section we estimated to be 
0.006 dex (see Table[2]i. In order to estimate how this error affects 
the measure of A A ge, we iterated the same procedure previously de¬ 
scribed, using for POPa the same set of isochrones, and for POPb 
two additional sets of isochrones with the same [Fe/H]=-0.67 and 
[o-/Fe]=+0.4 as for POPa, but with two different helium enhance¬ 
ments: AT = 0.029 + 0.006 = 0.035 and AT = 0.029 - 0.006 = 
0.023. We found that an error of 0.006 dex in AT translates into an 
uncertainty in Aage of ct AAG e(T) = 40 Myr. 

The AGB and FRMS models predict that the stars of 
SG have the same [Fe/H] and [o-/Fe] of stars of FG, as ob- 
served in many GCs. Mo reove r the spectroscopic results by 
IFeltzing. Primas & Johnsonl |2009|) confirm that in NGC 6352 the 
stars of the two populations should have the same metallicity 
(within 0.02 dex). In any case, we also considered the cases 
in which POPb has different [Fe/H] or [or/Fe] with respect to 
POPa, to study the impact of these variations on A A ge- As first 
test, we considered the case in which POPb has metallicity 
[F e/Hlpppb = [Fe/H|pop a ± 0.02 dex, consistent with the results 
by IFeltzing. Primas & Johnsonl J2009l l. We performed the proce¬ 
dure described above using for POPa a set of isochrones with 
[Fe/H]=-0.67, [ff/Fe]=0.4 and primordial helium, and for POPb 
[Fe/H]=-0.65, [tr/Fe]=0.4, and T = 0.285. We found that the 
relative age between the two populations is A A ge = 210 Myr. 
In the same way, we considered for POPb a set of isochrones 


Parameter 

Variation 

0"AAGE 

Photometry 

AT 

±0.006 

±110 Myr 
±40 Myr 

[Fe/H] 

±0.02 

±200 Myr 

[or/Fe] 

±0.02 

±150 Myr 

C+N+O 

<9AGE/<9[CNO] ~ 

-3.3 Gyrdex -1 


Table 3. Summary of all the uncertainties on Aage 


with the same parameters, except the metallicity that we set at 
[Fe/H]=-0.69. In this case we found Aage = -190 Myr Therefore, 
a variation on the metallicity of the second population of <5[Fe/H]= 
±0.02dex translates in a variation on A A ge of cr A AGE([Fe/H]) = 
±200 Myr. 

We then considered the impact of the ^-enhancement varia¬ 
tions on the estimate of A A ge- Similarly to what was done already 
for the variations in metallicity, we considered a set of isochrones 
for POPa with [Fe/H]=-0.67, [or/Fe]=0.4 and T = 0.256, and 
for POPb the same [Fe/H], enhanced in helium by AT = 0.029 
and d[o-/Fe] = ±0.02 with respect to POPa. We found that a 
variation of <5[o'/Fe]= ±0.02 dex leads to a variation on Aage of 
cr AA GE([a/Fe]) = ±150 Myr. 

It is possible that POPa and POPb have different C±N±0 
abundances. Unfortunately, we are not able to directly verify 
the eff ects of C±N ±Q abundances variations on the estimate of 
Aage - iMarino et al.1 d2012al ) presented a [Fe/H]-independent rela¬ 
tion between the C±N±0 variations and the age of a population 
<9AGE/<9[CNO] ~ -3.3Gyrdex If the POPb is CNO-enhanced 
with respect the POPa, then POPb is younger compared to the case 
in which both populations have the same CNO abundances, and 
therefore A A ge is larger. However we note that the m F8l4W and 
( m F606w - ™F8i4w) magnitude and color on which our analysis is 
based are the least affected by C, N, O variations. 

In Table [3] we provide a summary of all the uncertainties on 
Aage introduced by the variation of the considered parameters. 


7 CONCLUSIONS 

In this paper we have presented an analysis of the CMDs and two- 
color diagrams from WFC3/UVIS F275W, F336W and F438W, and 
ACS/WFC F625W, F658N, F606W and F814W photometry. 

We identified two stellar populations (named POPa and 
POPb), that, with appropriate combinations of magnitudes and col¬ 
ors, we could clearly distinguish on most of the evolutionary se¬ 
quences of the CMD, from the MS, to the SGB, RGB, AGB, and 
HB. 

Using a multi-color analysis method we have already applied 
to half a dozen of other GCs, we estimated a AT = 0.029 ± 0.006 
between the two populations. 

We also developed a new procedure for evaluating the dif¬ 
ference in age between the two populations hosted by NGC 6352. 
We considered the filter combinations that are the least affected by 
light-element variations. We studied the case in which POPa and 
POPb have the same [Fe/H], as i nferred from high-resolution spec¬ 
troscopy bv lFeltzing, Primas & Johnsonl l2009i ). but POPb has he¬ 
lium abundance enhanced by AT = 0.029 ± 0.006. We found that 
the two populations have a difference in age between the POPa and 
POPb of Aage = 10 ± 110 Myr, where the error includes the pho¬ 
tometric error and the error on the fit. The error in AT leads to an 
error on the relative age of <t A age = ±40 Myr, and a total error of 
120 Myr. This is the best relative ages between the two populations 
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we can estimate assuming that the two populations have the same 
[Fe/H] and [ff/Fe], 

In the case that POPb has a different [Fe/H] content 
of ±0.02 dex (as su g gested by high-resolution spectroscopy by 
IFeltzing. Primas & Johnsonll2009h . the relative ages would change 
by ±200 Myr. For a difference in [ff/Fe] of ± 0.02 dex, we found a 
A A ge variation of ±150 Myr. As shown bv lMarino et alj l l2012ah . 
if POPb is CNO-enhanced then it must be younger since the turn¬ 
off brightness is the same. Differences in metallicity between two 
populations hosted by a normal GC, however, are not predicted by 
“AGB” and “FRMS” models of formation of multiple stellar popu¬ 
lations and are also not observed for this cluster. 

In literature, there are numerous works on the relative age of 
multiple stellar populations hosted by anomalous GCs for which 
a spread in metallicity is measured. Surely , the most co ntroversial 
case is represented by u> Centauri. ISollima et alj ( 120081) identified 
four coeval pop ulations with different metallicities and helium con¬ 
tents, while lVillanova et"aD d2014l ) identified a large spread in age 
for each population and found that the most metal-rich population 
also is the oldest one, concluding th at at Centauri i s the re sult of 
a merger of two different progenitors. iMarino et aid d2012td ) found 
that the populations of M 22, characterize d by different me tallici- 
ties, are almost coeval within ~ 300 Myr. iLee et all d2013l ) found 
that the GC NGC 2419 hosts two populations characterized by 
a large differenc e in [ Fe/Hl, Y and with a relative age of about 
2 Gvr. lRoh et aid d201 ill compared models and CMDs of NGC 288 
to put constraints on the second stellar population hosted by this 
GC. They found that, in order to properly reproduce the observed 
CMDs, the second generation must be moderately metal enhanced 
by 0.16 dex, helium-enhanced by 0.03 dex, and younger by 1.5 Gyr 
than the first stellar generation. 

The work we present in this paper is different from previ¬ 
ous studies on the relative ages of multiple stellar populations in 
GCs, in the sense that this is the first attempt to measure relative 
ages within a multiple population GC with no observed signature 
of [M/H] dispersion. Assuming no difference in [Fe/H], [or/Fe], 
and C±N±0 content, the two populations of NGC6352 have a 1 cr 
spread in age ~ 120 Myr. A combination of small differences in 
[Fe/H] and [cr/Fe] of 0.02 dex would rise the total uncertainty on 
the relative age to ~280 Myr (within which the two populations are 
still coeval). 
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